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1 01 THE COSMIC SAT DIFFUSION IS A VIOLENT INTERSTELLAR MEDIUM 

A.M.Bykov and I.N. Toptygin 

M. I. Kalinin Poly technical Institute, Leningrad 195251* USSR 

1. Introduction. A variety of the available observational 
data on the oo anile ray (CR) spectrum, anisotropy and compo- 
sitlon are in good agreement with a suggestion on the diffu- 
sion propagation of CR with energy below 10(15) eW in the 
interstellar medium (see Ginsburg and Ptuskin, 1983)* The 
magnitude of the CR diffusion coefficient and its energy de- 
pendence are determined by ISM magnetic field spectra. 

Direct observational data on magnetic field spectra are 
still absent. In this work we present a theoretical model to 
the turbulence genezation in the multiphase ISM. Our model 
is based on the multiple generation of secondary shocks and 
concomitant large-scale rarefactions due to supernova shock 
interactions with interstellar clouds. We derive (also see 
Bykov and Toptygin, 1985) the distribution function for ISM 
shocks with account for supernova statistics, diffuse cloud 
distribution and various shock wave propagation regimes. 

This allows us to calculate the 191 magnetic field fluctua- 
tion spectrum and CR diffusion coefficient for the hot 
phase of ISM. 

2. The I9M turbulence model. Current observational data on 
the structure of I at turbulence have been summarized in the 
review by Armstrong et al.(1981). They have concluded that 
a variety of data obtained by different methods do not con- 
tradict to the assumption on the presence of continuous 
spectrum of ISM inhonogenuitles in a wide scale range. The 
energy density dependence on wavenumber k takes the form . 
W(k )<- k“<* , where l,4£oC$1.8 for 10(-20) 4 k < 10(-8) om** x 
(Armstrong et al., 1981). 

The ISM is characterized by a wide variety of phase star 
tea, so that the ISM turbulence is highly nonunifona. The 
possibility of developing an inertial range of weak MHD tur- 
bulence produced by energy cascading from long wavelengths 
to short wavelengths under the ISM conditions has been dis- 
cussed by Mclvor (1977). This discussion is based on the as- 
sumption about oollisional damping of MHD waves. We remark 
(Bykov and Toptygin, 1983) that the weak MHD turbulence 
spectrum can be developed by the energy cascading if only 
* the linear Landau damping of collisionless MHD waves is ef- 
fectively suppressed. The above suppression can be realized 
due to the distortion of the thermal plasma distribution 
function by finite amplitude waves ( 5"B£.0,1£ ) with seales 
1 3*10(17; cm. If such waves are absent, the® energy casca- 

ding of weak MHD turbulence seems to be highly unlikely as 
well as the appearance of some inertial range in the hot ISM. 

On the other hand, the fluctuations of ISM magnetic 
field and electron density may be connected with the presen- 
ce of strong turbulence. Here we treat the strong turbulence 
as an ensemble of noncoherent weakly interacting nonlinear 
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waves. For sueh a system the energy redistribution between 
wave-scales are determined by the evolution of nonlinear 
waves. The strong turbulence of Bangmuir solitons has been 
discussed by Budakov and fsytovleh (1978)., 

In the present paper we propose the model for a strong 
ISM turbulence. The main structural elements of our model 
are nearly spherical oollisionless shock waves with rarefac- 
tionsregions in the inner parts of the spheres. Sueh struc- 
tures may be formed due to multiple interactions of strong 
supernova shocks with interstellar clouds. The total energy 
input to the ISM turbulence in our model can be estimated 
as 10$ from the kinetic energy of large-scale motions pro- 
duced by superaovae and stellar wind of OB stars (about 
10(-27) erg* om(-3)* s(-l)). low we obtain the distribution 
function for ISM shocks of different strengths (also see 
Bykov and, Top ty gin, 1985). 

3. ISM shock wave distribution function. Theoretical consi- 
deration of an average distribution of ISM shocks is impor- 
tant to study the ISM structure (McKee and Ostriker, 1977 J 
Heathcote and Brand, 1983) and OR propagation (Blandford 
and Ostriker, 1980; Axford, 1981; Ginzburg and Ptuskin,1984). 
We define the shock strength distribution function P(m) as 
an average number of passings of shocks with Mach nuMoers 
IM through an arbitrary point of ISM disk per unit 
time (cf. Blandford and Ostriker, 1980). Supernova explosi- 
ons with local rate S (~10(-13) pc(-3)/year) produce strong 
primary shock waves in the ISM. Hie primary shock evolution 
is described by self similar spherical solutions with radius 
RQO - R„ 

where R 0 depends on w 

meters. The exponent is equal to <k* 2 for the Sedov so- 
lution; (a ■ 4,5 for the McKee and Ostriker solution, and 
Cn * 1,2 for the snow-plow model. 

It is known that an average number of passings of pri- 
mary shocks with Maoh numbers ^ M through an arbitrary 
point of ISM per unit time is given by (see e.g* Blandford 
and Ostriker, 1980; Arford, 1981) 

M *</*> ■ •’Hy*. 

Thus, P (M) » -dF/d /4 . Primary shock propagation through 
ISM with embedded diffuse clouds will be accompanied by ge- 
neration of secondary shocks (e.g. see SPitzer, 1982). The 
number of secondary shocks produo ed by a primary shock will 
be about 10(4), provided the cloud number density 

*5‘10(-4) pc(-3), and R 10(2,3) pc. Hie distribution fun- 
ction of secondary shocks (without account for secondary 
shock reflection from clouds) takes the forms 


OJLX OJi 

, ( 1 ) 
xhe energy release and ambient gas para- 


^5 

i r 



(*) 
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where is the diffuse cloud filling factor in the ISM. 
This formula takes account for shock reflection from dense 
sphere and the secondary shock propagation regime. For simp- 
licity, we consider a perfect gas with £ ■ 5/3 and take into 
account that secondary chock will be weak (i.e. 1<1) if 

the primary shock hasJA 4 fc 3 (for more detailed discussion see 
Bykov and Toptygin, 1985). 

As a result, we derive an approximate distribution func- 
tion with account for the weak secondary shock contribution* 






where C(3;4,5) 4,l-10(-3). It should be noted the strong 

dependence of P on M determined by the secondary shock con- 
tribution. Among the waves with JA 1 1,2, secondary shocks are 
dominant. Multiple interactions of secondary chocks with 
clouds can be taken into account by iterations of eq. (3)» 
i.e. by series expansion of P(JA) with respect to f 0 ^. 

4. ISM fluctuation spectra . In the discussed model nearly 
aphft-rift tti ahoctWw with rarefaction regions in the inner parts 
of spheres represent well-defined structural cells. She ISM 
velocity, density and magnetic field fluctuations are deter- 
mined by noncoherent overlapping of a large number of such 
structures. For statistically homogeneous scad isotropic en- 
semble of chocks (with strength jK ) and concomitant rarefac- 
tion the spectral energy density takes the form* 

uHk,^= h/fikj'f/M'fy kVfjwiT 1 (5) 

whichis valid for scales k$<><l, where S' is the shock front 
thickness. For oollisionless weak shocks in the hot ISM o’ £ 
10(12) cm. Bq. (5) quite satisfactorily extrapolates the well 
-known relations at kR»l and kR«l (Bykov and Toptygin, 1985). 

To evaluate the spectral energy density, W(k), one should 
average Eq. (5) over an ensemble of shocks with different 
strengths.^This yields W(k) * nr’(k) ♦ W (a <k), where 

W |l '(k) * 

Let us pVesent the obtained expression for W(k) in the short 
- wavelength limit (kR V>1) for the McKee and Ostriker (1977) 
model (d>3): 

w(k) = ( h si* c i* •, \ k 

R«(, being the mean radius of clouds; (^-1) * 10(-2) ; C(2;4»5) 
^ 9,4*10(-3). The long-wavelength spectral range has been 
considered earlier by the authors. The magnetic fluctuation 


70 


0G 7.2 - 14 


spectrum appears to be similar to the velocity fluctuation 
spectrum given by Eqs. (6), (7). 

Thus, Eqs* (6) determine the fluctuation spectrum in the 
model of turbulence produced by shocks from supernova explo- 
sions. The main mechanism of energy redistribution over wave- 
lengths is that due to the generation and subsequent self- 
similar evolution of secondary shocks. 

5. The OR diffusion . The magnitude of the CR diffusion coef- 
ficient along the field lines and its dependence on the 
energy E of energetic particles are determined by the magne- 
tic fluctuation spectra (e.g. Toptygin, 1985), 

*n * , (8) 
where the effect of large-scale magnetic fields on the par- 
ticle scattering at is taken into account. According 

to our result (7) for the McKee and Ostriker (1977) model, 
the spectral turbulence index V « 2. Hence, it follows from 
Eq. (8) that the CR diffusion coefficient along the regular 
magnetic field is energy-independent as long as E eR B • 

For the above mentioned parameters of the model we have 0 

pQ p ^ 

<X|i ~ $* 10 cm s x which is in satisfactory agreement with 
observational data on CR anisotropy and abundance. 
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